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Short title: INFLUENCE OF AIRMASS HISTORIES 



Abstract. The Goddard trajectory chemistry model was used with ER-2 aircraft data to test our 
current knowledge of radical photochemistry during the POLARIS (Polar Ozone Loss in the Arctic 
Region In Summer) campaign. The results of the trajectory chemistry model with and without 
trajectories are used to identify cases where steady state does not accurately describe the 
measurements. Over the entire mission, using trajectory chemistry reduces the variability in the 
modeled NO x comparisons to data by 25% with respect to the same model simulating steady state. 
Although the variability is reduced, NO x /NO y trajectory model results were found to be systematically 
low relative to the observations by 20-30% as seen in previous studies. Using new rate constants for 
reactions important in NOy partitioning improves the agreement of NO x /NO y with the observations but 
a 5-10% bias still exists. OH and H0 2 individually are underpredicted by 15% of the standard steady 
state model and worsen with the new rate constants. Trajectory chemistry model results of 0H/H0 2 
were systematically low by 10-20% but improve using the new rates constants because of the explicit 
dependence on NO. This suggests that our understanding of NO x is accurate to the 20% level and 
HO x chemistry is accurate to the 30% level in the lower stratosphere or better for the POLARIS 
regime. The behavior of the NO x and HO x comparisons to data using steady state versus trajectory 
chemistry and with updated rate coefficients is discussed in terms of known chemical mechanisms and 


lifetimes. 
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Introduction 

Ozone abundances outside the tropical latitudes show a distinct seasonal cycle. In the northern 
high latitude stratosphere, column ozone abundances decrease from a maximum in late spring to a 
minimum in early fall [Newman et al., 1997], To predict trends in ozone and responses of the 
atmosphere to perturbations such as emissions of a proposed fleet of supersonic transport aircraft 
[Stolarski et al., 1995], we must understand the processes driving the high latitude seasonal cycle in 
ozone concentrations. 

During the summer months, transport of high ozone mixing ratios from low latitudes and high 
altitudes to the high latitude lower stratosphere is small [Wu, 1987], Ozone production by photolysis 
of O z is also small relative to tropical latitudes [Johnston, 1975], With summertime stratospheric 
temperatures being well above PSC thresholds [Rosenlof, 1996], ozone loss due to chlorine activation 
on polar stratospheric clouds is unlikely. 

The seasonal loss in the northern high latitudes between late spring and fall is generally 
attributed to in situ photochemical reactions [Perliski, 1989]. The main chemical loss cycles are due to 
NO, (NO and NO,), HO, (OH and H0 2 ) with minor contributions from the chlorine and bromine 
catalytic cycles (in conditions where aerosol is not enhanced by recent volcanic activity). Catalytic 
loss cycles begin with reactions of NO, OH and halogen atoms with ozone, followed by the loss of an 
additional odd oxygen (0 3 and O) and reformation of NO, OH, Cl, or Br: 

NO + 0 3 -> N0 2 + 0 3 (Rl) 

NO, +Q -> NO + O, (R2) 

0 3 + O -> 0 2 + 0 2 

OH + 0 3 -> HO, + O, (R3) 

HO, +0, -> OH -f 20, (R4) 

0 3 + 0 3 -> 30 2 



(R5) 

(R6) 
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X + 0 3 -> XO + O, (R5) 

XO + O -> X + a (R6) 

o + o 3 -> o 2 + 0 2 . 

where X = Cl or Br. The HO x cycle is coupled to NO x and halogens by the following reactions: 


H0 2 + NO -> N0 2 + OH (R7) 

H0 2 + XO -> HOX + 0 2 (R8a, b) 

HOX + hv -> OH + X (R9). 


HO x partitioning also depends on CO 

0H + C0->C0 2 + H (RIO) 

H + 0 2 + M->H0 2 + M (Rll). 

In situ photochemical loss in summertime is primarily due to the NO x catalytic cycle [Perliski, 1989], 
Figure 1 shows the interrelationship of the NO x species, NO y reservoirs, and reaction pathways 
between them. During the spring and summer, increased solar illumination prevents the formation of 
N 2 0 5 from N0 2 + N0 3 by rapidly photolyzing N0 3 , resulting in higher abundances of NO and NO> 
Figure 1 also shows that species from the other catalytic cycles such as OH and CIO are interrelated 
with the NO x cycle. Understanding the reactions and pathways of the NO x cycle is important not only 
for assessing ozone loss due to NO x but also the influence of NO x on the HO x and C10 x catalytic loss 
cycles. 


The POLARIS campaign provides the first ER-2 opportunity to test the photochemistry 
associated with high latitude ozone loss between spring and fall. At high latitudes in the summer, 
airmasses experience continuous or near continuous sunlight. Under these conditions, radical species 
with short lifetimes such as NO x , HO x and C10 x are often assumed to be in steady state. Steady state 
model calculations have been used in previous ER-2 campaigns to model the diurnal variation of 
radical species [Salawitch et al., 1994a, b] during the 1993 NASA Stratospheric Photochemistry, 
Aerosols, and Dynamics Expedition (SPADE) and in the investigation of NO y partitioning [Gao et al., 
1997] during 1994 NASA Airborne Southern Hemisphere Ozone Experiment/Measurements for 
Assessing the Effects of Stratospheric Aircraft (ASHOE/MAESA). However, the data in some cases is 
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not accurately described by steady state conditions. Kawa et al. [1993] found that for one high- 
latitude case, a photochemical model that accounts for latitude and temperature changes 10 days prior 
to sampling was needed to provide better agreement with observations. In this paper we show how 
the Goddard trajectory chemistry model can be used to identify POLARIS special cases, put them on a 
consistent basis with the steady state cases, and in doing so help to separate out systematic biases 
linked to model chemical mechanisms or rates. Our findings show that use of the trajectory chemistry 
model reduces the variability of calculated NO x relative to the steady state model. As a result of the 
interrelationship between NO x and HO x , HO x variability is reduced also. Trajectory chemistry model 
results of NO x are found to be systematically low with respect to the observations using rates from the 
JPL-97 compendium, consistent with the deficit seen in Gao et al. [1999]. Because the NO x and HO x 
cycles are linked, the NO x bias must be accounted for before systematic offsets in the HO x model 
results can be assessed. Use of the new rate constants of Brown et al. [1999a, b] improves modeled 
NO x with observations but agreement of OH and HO, individually with observations worsens slightly. 

Model Description 

The chemistry on trajectory model consists of two parts, a trajectory module [Schoeberl et al., 
1993] and a chemistry module [Kawa, et al., 1997], Using the Goddard assimilated meteorological 
fields [Algorithm and Theoretical Basis Documents of the Data Assimilation Office, 1997], air parcels 
intersecting the ER-2 flight path every 432 seconds (-86 km) are advected backwards 10 days. 

Latitude, longitude, temperature, pressure and solar zenith angle are calculated at 15 minute time steps 
along the isentropic back trajectory. 

The chemistry portion of the model includes 30 species and 120 reactions with reaction rates 
from the JPL-97-4 compendium [DeMore et al., 1997]. The chemical mechanism is designed for 
stratospheric conditions. Constituent mixing ratios of long-lived species NO r H,0, CH^, 0 3 and CO 
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are initialized from ER-2 measurements [Gao et al., 1997; Hintsa et al., 1998; Webster et al. ( 1994; 
Proffitt et al., 1989]. Cl y and Br y were estimated from measured N 2 0 and CFC1 3 and the correlations 
of Woodbridge et al. [1995] and Wamsley et al. [1998], Initial Cl y partitioning is set by measured 
HCI with the balance in C10N0 2 . Initial CIO concentrations are set to zero. Initial NO y partitioning 
is set by measured NO x , C10N0 2 from Cl y , small fractions in N 2 O s (0.5% of the N0y-C10N02-N0x) 
and H0 2 N0 2 (3% of NO y -ClONO 2 -NO x -0.1*Br y -N 2 O 5 ) and the remainder in HN0 3 . Initial CH 3 OOH is 
from the Goddard 3D CTM interpolated to the trajectory starting point’s latitude and longitude 
coordinates and H 2 0 2 is determined from a fit to measured 0 3 . Other species are initially zero or 
partitioned in instantaneous steady state. 

Heterogeneous chemistry parameters are established assuming measured particles are composed 
of a ternary solution of sulfate, nitrate and water [Kawa et al., 1997], The parameters are used 
primarily for the N 2 0 5 and BrON0 2 hydrolysis reactions whose reaction sticking coefficients are 0.1 
and 0.4, respectively, for this analysis. For POLARIS conditions the HN0 3 content of the particles is 
negligible. H 2 S0 4 mixing ratio is obtained from the Focused Cavity Aerosol Spectrometer FCAS 
[Jonsson et al., 1995] volume by assuming vapor pressure equilibrium [Carslaw et al., 1995], HN0 3 
varies with the chemistry along the back trajectory while H 2 S0 4 , measured aerosol number, and H 2 0 
mixing ratios are assumed to be conserved along the back trajectory. Temperature variations along the 
trajectory produce composition changes leading to a change in the total particle volume. Reactive 
surface area is calculated from measured particle number and calculated particle volume. 

Photolysis rates are obtained from a table look-up based on albedo, overhead ozone, solar 
zenith angle, pressure and temperature. The tabulated calculations are from the radiative transfer 
model of Anderson and Lloyd [1990] and Anderson et al. [1995]. Albedo is determined along the 
back trajectory by interpolating total ozone mapping spectrometer (TOMS) reflectivities to the 
trajectory’s latitude /longitude coordinates. Overhead ozone is determined by scaling climatological 



profiles and TOMS data along the trajectories to match overhead ozone from the Composition and 
Photodissociative Flux Measurement (CPFM) instrument [McElroy, 1995] at the flight track. 

The model can be also run fixed in latitude and longitude at the point where the trajectory 
intersects the flight path. Models run at a fixed point are initialized in the same manner using flight 
data and the diurnal chemistry proceeds at that point for 20 days. Running the model in this mode is 
considered to approach steady state as constituents change less than 5% over the last ten days of the 
model run. We designate the different modes as trajectory chemistry and point chemistry in the 
discussion of results. 

Comparison with PSS model 

The photochemical stationary state (PSS) model used in previous ER-2 aircraft campaigns and 
in POLARIS assumes that each species reaches a balance between production and loss over 24 hours 
for the temperature, pressure and latitude of a selected point along the flight track integrated over a 
diurnal cycle [Salawitch, 1994a]. Measured NO y , 0 3 , H 2 0, and CH 4 are used to initialize the model, 
and Cl y and Br y are inferred from relationships with CFC and brominated source gases [Woodbridge et 
al., 1995; Schauffler et al., 1993]. There are 35 reactive species and 200 chemical reactions with rate 
constants and cross sections based on the JPL-97 compendium. Photolysis rates are calculated using a 
radiative transfer model that accounts for Rayleigh and aerosol scattering [Prather, 1981; Salawitch et 
al., 1994a]. Total column ozone and ultraviolet albedo are specified from observations of TOMS or 
CPFM measurements when available. 

Figure 2a is a timeseries plot of NO, N0 2 and NO x /NO y observations and model results using 
JPL-97 rates from the point chemistry model and PSS for the POLARIS flight on July 7, 1997. 

Steady state conditions were dominant since solar zenith angles were between 55 and 70° and back 
trajectories show air parcels experience continual sunlight, indicative of high solar exposure and short 
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N0 X lifetimes (< 1-2 days). The timeseries of NO, N0 2 and NO x /NO y shows that the trajectory 
chemistry model run in the point chemistry mode agrees well with the PSS model results. Point 
chemistry model results of NO x agree within 5% of the PSS results. The 5% underestimate of NO x in 
the point chemistry model is due to the underestimate of N0 2 . These small differences are traceable 
to small differences in the albedo and photolysis rates but agreement between both models shows that 
both models calculate a similar photochemical steady state over a wide latitude range (5-90°N). 

Figure 2b is a timeseries plot of OH, H0 2 and 0H/H0 2 observations and model results similar 
to figure 2 for July 7, 1997. The First two panels show that HO x in the point chemistry model agrees 
with the PSS model results to within 15%. Although not shown here, this HO x point chemistry 
agreement with PSS mode results is consistent for all the POLARIS. HO x overall, however, is 
underestimated relative to the observations as seen in the third panel. It is important to recognize that 
because of the interrelationship between NO x and HO x cycles, the bias seen in NO x influences HO x and 
is addressed later. 

We extend the analysis between point chemistry and PSS to include flights over the entire 
mission. Figure 3a shows the scatter plot comparison of modeled NO x /NO y using JPL-97 rates from 
the point chemistry model and PSS. The flights of 970630 and 970710 were not included since they 
had very few trajectory points above 90 hPa to analyze (< 5). The NO x /NO y ratio in the point 
chemistry model compares to within 3% of the PSS model with a slope and correlation coefficient 
near 1.0. The point chemistry model can now be used as a reference to steady state to begin 
examining cases where point chemistry/steady state does not accurately describe observations during 
POLARIS by comparing results with the trajectory model results. 

Figure 3b compares modeled NO x /NO y using trajectory chemistry to point chemistry. In 
general, trajectory chemistry model results agree with or lie above the 1:1 line. However, portions of 
the flights of April 24 (blue diamonds), April 26 (green triangles), September 8 (green pluses) and 
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September 21, 1997 (yellow triangles) lie more than 20% off of the 1:1 line. These flights represent 
cases where point chemistry does not accurately describe the NO x observations. 

NO x Example flight 

The flight of April 26, 1997 during the spring deployment presents an interesting case study 
for use of trajectory chemistry. The ER-2 flew out of Fairbanks (65°N) to the pole and back with a 
profile at the pole. Solar zenith angles during the flight were between 55 and 75 degrees and the 
aircraft sampled the edge of the persistent winter vortex as seen in reduced amounts of CH 4 , N 2 0 and 
CO [Webster et al., 1994]. Figure 4 shows the NO, N0 2 and NO x /NO y timeseries for the trajectory 
and point chemistry model results using JPL-97 rates and observations. The latitudes sampled during 
the flight are shown on the top axis. At first glance, it appears that the observations are in steady state 
as all three panels show that the point chemistry model is in reasonable agreement with the 
observations between 64000 and 80000 UT seconds. The trajectory chemistry model does not agree 
well with the observations, underestimating NO x . Figure 5 shows 10-day back trajectories of the 
parcels intersecting the flights path’s on the outbound leg of the flight to the pole. Air parcels 
experienced a strong southerly advection resulting in a large latitudinal excursion over the last few 
days prior to ER-2 sampling. In such a situation, it is quite likely that the parcels are not in 
photochemical steady state. 

Figure 6 shows the solar zenith angle, NO x , and HN0 3 histories in the point and trajectory 
chemistry models for the ten days prior to sampling of parcel 26, marked with an arrow in Figure 4 as 
P26. In the trajectory model, the air parcel spends time in darkness (SZA > 95) while in the point 
chemistry model it does not experience darkness at all. The trajectory chemistry model accounts for 
latitude variations along the back trajectory while the point chemistry model holds the point fixed with 
no latitudinal variation. The longer time in darkness along the back trajectory results in less photolysis 
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of HN0 3 , more formation of N 2 0 5 , and consequently less NO x . 

For the flight of April 26, 1997, the steady state/point chemistry model agrees well with 
measured NO x but for the wrong reason. The difference between the trajectory chemistry model 
results with observations is a better indicator of the systematic difference between model and actual 
photochemistry than the point chemistry model, even though the point chemistry model is in closer 
agreement . Air parcels experienced changes in the solar exposure on time scales much shorter than 
the time scale for the photochemistry to relax to steady state. Previous steady state modeling results 
by Gao et al. [1997, 1999] have shown the PSS model to underpredict NO x /NO y by 20-30% similar to 
the comparison with the trajectory chemistry in Figure 4. Comparison of NO/NO y from the chemistry 
model with all of the POLARIS observations (below) shows that NO x /NO y is consistently 
underestimated by a similar amount. 

Variability studies 

The trajectory chemistry model complements and improves use of the steady state assumption 
in modeling aircraft data by identifying flights where trajectories influence photochemistry. In these 
flights, the steady state assumption may not be adequate in modeling those regions. In figure 7, we 
plot the comparison of the trajectory and point chemistry model results using JPL-97 rates versus all 
mission observations for NO x /NO y . The solid black line represents the one-to-one line and the solid 
red line denotes a linear least squares fit to the data. The trajectory and point chemistry model results 
display similar 20-30% deficiency in NO x relative to the observations indicating that there is a 
systematic bias associated with the model photochemistry. A closer look at the plots shows that the 
variability of the trajectory chemistry NO x /NO y results about the fit line is less than the variability of 
the point chemistry model results by -30% as seen in the differences in the standard deviations about 
the fit line (0.019 vs 0.013). Since NO x has a photochemical lifetime of - 1 week, it is still adjusting 
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to changes in photolysis caused by changing latitude. This source of variance in the model-data 
comparisons is removed using the trajectory chemistry. The analysis shows that, while POLARIS 
NO x /NO y data as a whole may be reasonably well-characterized by steady state, individual cases can 
vary by 25% or more. By reducing the variability, we are able to better diagnose systematic chemical 
relationships and assess our understanding of the photochemical processes associated with ozone loss. 

Model calculations with new rates 

We have seen that trajectory chemistry reduces the NO x variability but exhibits a systematic 
bias in the photochemistry. During POLARIS, the increasing solar illumination and warmer 
temperatures prevent N 2 0 5 formation and N 2 0 5 hydrolysis from controlling the NO x /NO y partitioning 
[Farman et al., 1985]. Under these conditions, the NO x /NOy ratio is controlled primarily by the 
following reactions: 

NO, + OH + M -> HNO3 + M (R1 1) 

OH + HN0 3 -> N0 3 + H 2 0 (R12) 

HN0 3 + hv -> OH 4- N0 2 (R13). 

Recently, Brown et al. [1999a, b] reported rates for R1 1 that are 10-20% slower than recommended in 
JPL-97 and rates for R12 that are faster by 40-100% for stratospheric temperatures and pressures. 
Figure 8 shows the result of incorporating these rates into the trajectory chemistry model for the April 
26, 1997 flight. Both NO and N0 2 are increased by approximately 15-25% and are in better 
agreement than the model using JPL-97 rates for reactions R11-R13 above. Figure 9 plots the 
comparison of modeled NO x /NO y using Brown et al. and observed NO x /NO y for the entire mission. 
The solid black line through the data is the least squares fit. Comparing Figure 9 and the top panel of 
Figure 7, inclusion of Brown et al. [1999a, b] increases NOx by 25% and reduces the bias between 
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modeled NO x /NO y and observations by ~70% . More NO x is available as a result of slower loss of 
N0 2 from reaction Rll and faster production from R12. The agreement of modeled NO x /NO v in 
figure 9 is similar to the improved agreement seen in Gao et al. [1999]. 

Inclusion of the updated rate constants improves the agreement with observations but a 5-10% 
bias still exists. Instrumental uncertainties for NO, N0 2 and NO y [Gao et al., 1997] collectively are on 
the order of ±40% and so the bias is within experimental error. Other sources of error which may 
contribute to the bias include input albedo and 0 3 column above the ER-2 which are key parameters in 
establishing the model photochemistry. Albedo influences the partitioning of NO x between NO and 
N0 2 through J N02 (refer to Figure 1) but sensitivity studies show that doubling the albedo results only 
in a 10% increase in model NO x . Based on the analysis of TOMS irradiances by Herman and Celarier 
[1997], the albedo would not be expected to be uncertain by a factor of two. CPFM measurements of 
overhead ozone have an uncertainty of ±3% [McElroy, 1995] and trajectory model sensitivity studies 
for POLARIS conditions show a 3% increase/decrease in the overhead column results in a 3% 
decrease/increase in NO x abundance. An evaluation of the NOx production and loss sensitivity to 
changes in chemical concentration, rate coefficients and radiative conditions sensitivities during 
POLARIS by Perkins et al. [1999] demonstrated that changing the overhead column by 10% or 
decreasing the OH + N0 2 rate constant another 10% would be necessary to reduce the NO x deficiency. 
It is pointed out, however, that changes to NO x and NO y concentration measurements or photolysis 
rates may also contribute, thus making it difficult to separate out a single source. 

Effects on HOx 

The relationship between the NO x catalytic cycle and the HO x species implies that changes in 
NO x with the use of the new rates constants should lead to changes in the HO x . The HO x production 
and loss processes are best represented by OH since the majority of the HO x production and loss 
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processes involve the OH radical. HO x production through OH formation is primarily characterized by 
HN0 3 photolysis and reaction of H 2 0 with 0(‘D) followed by smaller contributions from HN0 4 and 
CH 2 0 photolysis at lower solar zenith angles [Hanisco et al., 1999] and contributions from HOBr 
photolysis at higher solar zenith angles [Wennberg et al., 1999]. The reactions of OH with NO^ OH 
with HN0 3 and reaction of H0 2 with N0 2 represent the primary HO x loss processes [Wennberg et al., 
1990, Hanisco et al., 1999]. 

Use of the Brown et al. rates will affect HO x directly by the decrease in the k 0H+N02 and 
increase in the ko H+ HN03 react i° n rates - The reaction of OH + N0 2 is roughly two times as significant 
as OH + HN0 3 in the loss of HO x . Since k 0H+HN03 increases by twice the amount that k 0H+N02 
decreases, the rates roughly cancel. HO x loss rates are affected indirectly by the increased N0 2 and 
decreased HN0 3 that result from using the new rate constants. The fractional change in N0 2 is large 
and tends to dominate over the small fractional decrease in HN0 3 so that the net effect is an increase 
in the HO x loss rates. Since the OH + N0 2 reaction makes up about half the total loss rate, a 20% 
increase in N0 2 using the Brown rate constants results in a 10% decrease in OH relative to JPL-97 
rates. H0 2 concentrations respond to the changes in the HO x loss rates but also to changes in the 
0H/H0 2 partitioning. Increased HO x loss rates result in lower H0 2 because there is less OH to 
convert to H0 2 and increased NO x results in lower H0 2 because H0 2 is converted into OH more 
quickly. 

The effects of trajectory chemistry and the use of the updated rate constants of Brown et al. 
[1999a, b] on OH, H0 2 and 0H/H0 2 can be seen in the statistical results for three model runs 
presented in Table 1. With the recommended JPL-97 rates, trajectory chemistry reduces the variability 
of OH by about 20% from 0.051 to 0.040 in the standard deviation about the fit as a result of the 
reduction in the variability of HN0 3 and N0 2 (see Figure 7). Using trajectory chemistry with the 
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Brown rate constants, the modeled OH is in lesser agreement than using the JPL-97 rates (-0.222 vs. 
-0.300) despite the fact that NO x /NO y is in better agreement with the observations (see Figure 9). This 
suggests that there is still a missing source or overestimation of a sink in the model’s HO x production 
and loss, consistent with Wennberg et al. [1999] who demonstrated underestimates of HO x but at 
higher solar zenith angles. Regardless of solar zenith angle, there is an underestimate of HO x . 

OH and H0 2 interconvert more rapidly than HO x production and loss processes so that OH and 
H0 2 are in photochemical steady state [Wennberg, 1990], Cohen et al. [1994] demonstrated that the 
OH/H0 2 ratio is well defined by the steady state equation 

[OH]/[H0 2 ] = k7fNQl + k4fO,l + k8alC101 + k8blBr01 (1) 

k3[OJ + klO[CO] 

where k3, k4, k7, k8, klO are the rate constants from reactions R3, R4, R7, R8 and R10. For the high 
latitude summer stratosphere, Lanzendorf et al. [1999] show that equation 1 simplifies to 
[OH]/[H0 2 ] = k7rNQl + k4rO,l (2) 

k3[0 3 ] 

Using the rates of Brown et al., the difference between the modeled and measured OH/H0 2 decreases 
(-0.211 to -0.139). Since the updated rate constants increase NO x , the 0H/H0 2 ratio increases also 
since it is governed only by the abundance of NO as ozone remains constant to first order. 

The effects of trajectory chemistry and inclusion of the Brown rate constants on H0 2 are more 
dramatic than those seen with OH. Using the JPL-97 rates, the average value of modeled HO, is very 
close to the observed HO, and not very different in going from point chemistry to trajectory chemistry 
(0.005 vs -0.006). However, trajectory chemistry reduces the modeled HO, variability more than it 
reduces the OH variability (0.61-0.41/0.61 33% vs 0.051-0.040/0.51 22%). Trajectory chemistry has 
more of an effect since H0 2 abundances are primarily controlled by NO x through the rapid 
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repartitioning reaction of H0 2 + NO -> OH + NO : . As the NO x variability improves, the H0 2 
variability improves as well. Replacing the JPL-97 rate constants with the updated rate constants of 
Brown et al. worsens the agreement between modeled and observed H0 2 (-0.006 to -0.185). The 
change in the modeled H0 2 with the trajectory chemistry using the Brown rate constants is greater 
than the change in the modeled OH (-0.222 to -0.300 vs -0.006 to -0.185) since the Brown chemistry 
affects HO x sources and sources as well as the partitioning between OH and H0 2 . 

The sensitivity of HO x to changes in model parameters is different than the sensitivity of NO x . 
While albedo is important in partitioning NO x into NO and N0 2 as shown by Schwartz et al [1999], 
sensitivity studies show that tripling the albedo from 0.2 to 0.6 only results in a 15% decrease in the 
HO x species. As discussed earlier, uncertainties in the albedo are unlikely to be of this magnitude and 
should not contribute to the modeled to observed HO x discrepancy. Sensitivity studies of CPFM 
measurements of overhead ozone show that a 10-15% increase/decrease in the overhead column results 
in a 10-15% decrease/increase in HO x since solar flux is important in the key HO x production 
processes involving the formation of 0( J D) from ozone [Schwartz et al., 1999] for reaction with H 2 0 
and CH 4 and photolysis of HN0 3 , and HN0 4 . It is important to note that due to the coupling between 
the HO x and NO x families, potential errors in HO x can affect NO x production and loss as well. NO x 
production is primarily characterized by HN0 3 photolysis and reaction with OH and NO x loss by the 
reaction of N0 2 with OH and BrON0 2 hydrolysis. Sensitivity studies by Perkins et al. [1999] have 
shown that a 10% increase/decrease in the OH concentration decreases/increases the production to loss 
ratio by 5%. Overall, modeled HO x agreement with observations is within the measurement 
uncertainties using the Brown rate constants. 



16 


Conclusion 

The POLARIS mission was designed to study chemical ozone loss processes which are 
accelerated as a result of increased solar illumination at summer high latitudes. NO x is primarily 
responsible for this stratospheric loss from late spring to early fall at high latitudes. Two useful tools 
for analyzing the radical abundances observed during POLARIS are the steady state and trajectory 
chemistry models. The trajectory chemistry model, initialized from observations and run in the fixed 
point mode agrees well with a steady state model, establishing a consistent basis from which we can 
begin to separate out the effects from trajectories and systematic biases in the photochemistry. An 
investigation of all mission flights revealed that as a whole POLARIS data could be adequately 
characterized by steady state but that trajectory history influences radical abundances significantly in 
several flights. The point chemistry model produces too much NO x as a result of longer, more 
consistent sun exposure than the shorter, more variable sun exposures encountered when trajectory 
history is taken into account. A comparison of the trajectory and point chemistry model results with 
the observations demonstrates that using the trajectory model reduces the variability in the NO x model 
results and is valuable in constraining relationships among the radical species. However, the trajectory 
chemistry model results consistently show a NO x deficiency representing a systematic bias in the 
photochemistry as in other steady state model results. The variability of HO x (OH and HO^ is 
reduced with trajectory chemistry because of the reduced variability in NO x and HN0 3 but OH is 
systematically low compared to observations by 20%. The NO x bias must be accounted for to 
properly assess the HO x chemistry in the model. 

To address the systematic bias in NO x , we included the new rate constants of Brown et al. 
[1999a, b] for OH + N0 2 and OH + HN0 3 in the trajectory chemistry model. Use of these rate 
constants improves the agreement of modeled NO x with observations to within the uncertainties of the 
measurements (± 40%) but a small bias (5-10%) still exists. The 0H/H0 2 ratio agreement with 
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observations improves with the new rate constants since the ratio is explicitly dependent on NO. 
Inclusion of Brown et al. rates increases the model offset of OH and HO, individually relative to 
observations, consistent with the underestimate in HO x seen by Wennberg et al [1999] and linked to a 
missing source in the current photochemistry used in models. The offsets of OH, H0 2 and 0H/H0 2 , 
however, are as a whole still within the 2-sigma measurement uncertainties of 25, 30, and 14%, 
respectively. 
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Figure Captions 

Figure 1. Schematic of the main reaction pathways for NO y species in the lower stratosphere. 

Figure 2. (a) Timeseries plot of NO, N0 2 and NO x /NO y for the flight of July 7, 1997. (b) Timeseries 
plot of OH, H0 2 , 0H/H0 2 for the flight of July 7, 1997. Latitude covered during flight is noted on 
top y-axis. Solid line denotes observations, dash-dot line represents the model results from the 
photochemical steady state model of Salawitch (1994a) using JPL-97 rates and the dotted line 
represents the model results from the Goddard trajectory model run in the fixed point mode using JPL- 
97 rates. 

Figure 3. (a) Scatter plot of GSFC point chemistry NO x /NO y model results versus JPL photochemical 
steady state model results for all mission flights. Data below 90 hPa excluded. Specific flights 
denoted in legend with the date in year, month day format (YYMMDD). Solid black line represents 
1:1 line and solid red line denotes linear least squares fit to the data, (b) Scatter plot of trajectory 
chemistry versus point chemistry NO x /NO y using JPL-97 for all mission flights. Solid black line 
represents 1:1 line and solid red line denotes linear least squares fit to the data. 

Figure 4. Timeseries plot of NO, N0 2 and NO x /NO y for the flight of 970426. Latitude covered 
during flight is noted on top y-axis. Small black dots denote observations, dotted line represents the 
model results from the Goddard trajectory chemistry model run at a fixed point (point chemistry) using 
JPL-97 rates and the dash-dot line represents the model results from the Goddard trajectory model 
using JPL-97 rates. The position in the time series of parcel 26 is denoted by the arrow P26 at 70500 


UT sec. 
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Figure 5. Back trajectory plot for points along the ER-2 flight track on the outbound leg from 
Fairbanks, AK (parcel 0) to the pole (parcel 30). Arrow heads mark trajectory positions every 24 
hours at OOz. 

Figure 6. Comparison between the trajectory and point chemistry model history of (a) solar zenith 
angle, (b) NO y species and (c) HN0 3 for parcel 26 marked in Figure 4. 

Figure 7 . Scatter plot of modeled NO x /NO y using trajectory chemistry with JPL-97 rates and point 
chemistry with JPL-97 rates versus measured NO x /NO y for the POLARIS mission. The legend of 
Figure 3 describes the figures for each flight. Solid black line represents the 1:1 line and solid red 
line represents the linear least squares fit to the data. 

Figure 8. Timeseries plot of NO, N0 2 and NO x /NO y for the flight of 970426. Latitude covered 
during flight is noted on top y-axis. Small black dots denote observations, dash-dot line represents the 
model results from the Goddard trajectory model using JPL-97 rates and the dotted line represents the 
results from the Goddard trajectory chemistry model with the rates of Brown et al. (1999a, b). 

Figure 9. Scatter plot of modeled NO x /NO y using trajectory chemistry with Brown et al. rates (1999a, 
b) versus measured NO x /NO y for the POLARIS mission. Solid black line through the entire plot 
represents the 1:1 line and solid line through the data represents the linear least squares fit to the data. 



Table 1. Statistics of modeled OH, H0 2 , and OH/HO, using point and trajectory chemistry with 
JPL-97 and Brown 1999 rates versus observed OH, H0 2 and OH/HO 2 . 


Point Traj Traj 

JPL-97 JPL-97 Brown 99 


OH 


Avg % diff (Mod-obs/obs) 

-0.220 

-0.222 

-0.300 

Slope 

0.842 

0.831 

0.762 

Y-int 

0.024 

0.019 

0.025 

Sdev about fit 

0.051 

0.040 

0.040 

Sdev about 1:1 

0.135 

0.134 

0.182 

H02 




Avg % diff (Mod-obs/obs) 

0.005 

-0.006 

-0.185 

Slope 

0.980 

0.978 

0.852 

Y-int 

0.097 

0.063 

0.087 

Sdev about fit 

0.612 

0.417 

0.348 

Sdev about 1 : 1 

0.615 

0.420 

0.716 


0H/H02 


Avg % diff (Mod-obs/obs) 

Slope 

Y-int 

Sdev about fit 


-0.213 

0.703 

0.013 

0.014 

0.046 


- 0.211 

0.679 

0.017 

0.013 

0.046 


-0.139 

0.824 

0.006 

0.016 

0.032 


Sdev about 1 : 1 
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10-Day Back Trajectories for 970426 
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